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Beschreibung 



Verfahren zur Bildbearbeitung ausgehend von einem unter Ver- 
wendung eines Kontrastmittels auf genommenen Computertomogra- 
5 phie (CT) -Bildes einer Lunge und CT-Gerat zur Durchfiihrung ei- 
nes solchen Verfahrens 

Die Erfindung betrifft ein Verfahren zur Bildbearbeitung aus- 
gehend von einem unter Verwendung eines Kontrastmittels ge- 
10 wonnenen Computertomographie (CT) -Bildes einer Lunge als Ori- 
ginalbild. Die Erfindung betrifft auSerdem ein CT-Gerat zur 
Durchfiihrung eines solchen Verfahrens . 

Die Computertomographie gewinnt fiir die Diagnose von pulmona- 
15 len Embolien (PE) zunehmend an Bedeutung. Zur Befundung wer- 
den dabei CT-Angiograpien der thorakalen GefaSe betrachtet. 

Zur Diagnose von PE wird eine CT Angiographie der Lunge 
durchgef iihrt . Dem Patienten wird Kontrastmittel mittels einer 

20 Kontrastmittelpumpe injiziert. Nach einigen Sekunden Pause 

wird ein Spiral-CT der Lunge durchgef iihrt , aus dem ein Stapel 
von in Richtung der Patientenlangsachse auf einanderf olgende, 
vorzugsweise aneinander angrenzende Korperschichten darstel- 
lenden axialen Bildern, d-h. Transversalschnittbildern, ge- 
wonnen wird. Die kontrastmittelgef tillten GefaSe sind auf CT- 
Bildern durch die Dichteanhebung deutlich zu erkennen. Die 
Befundung erfolgt an den axialen Bildern. GefaSe sind abhan- 
gig von der Orientierung der GefaSe zur Schnittebene als hel- 
le Linien oder helle Punkte im sogenannten Lungenf enster 

30 (Zentral-Wert des Fensters= -600 HU, Breite des Fensters = 

1500 HU) dargestellt. Beim „Blattern'* durch den Bilderstapel 
wird der gesamte GefaSbaum verfolgt und auf Verschllisse 
(Thromben) iiberpriift. Verschllisse sind als dunklere Materie 
im GefaS sichtbar. Nach dem Verschluss ist das GefaS nicht 

35 mehr, oder nur noch teilweise mit Kontrastmittel gefiillt und 
wird dadurch dunkler dargestellt. Folgende Probleme bestehen: 
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. - Verschliisse konnen libersehen warden, 

der hemodynamische Effekt der Thrombose kann nicht festge- 
stellt werden, 

5 

- Verschliisse an kleinen GefaSen (subsegmental) konnen im CT 
Bild durch die begrenzte Auflosung nicht sicher erkannt 
werden , 

10 - der Stenosegrad, d.h. die Schwere des Verschlusses kann 
nicht sicher bestimmt werden, und 




andere pathologische Veranderungen konnen irrtiimlich als 



thromboembolischen Verschluss diagnostiziert werden (z.B. 
15 verschlossene Bronchien) . 

Der Erfindung liegt die Aufgabe zugrunde eine zusatzliche 
Auswertemoglichkei t fur die bei der Diagnose von PE mittels 
Computertomographie anfallenden CT-Bilder anzugeben, die die 
20 Voraussetzung fiir eine erhohte Diagnosesicherheit bei PE bie- 
tet. Der Erfindung liegt auSerdem die Aufgabe zugrunde ein 
fiir diese Zwecke geeignetes CT-Gerat anzugeben. 

Nach der Erfindung wird diese Aufgabe gelost durch ein Ver- 
fahren mit den Merkmalen des Patentanspruchs 1. 

Infolge der Ermittlxing der Lungenparenchympixel und der Ein- 
farbung der Lungenparenchympixel, d.h. derjenigen Pixel des 
Originalbildes, die Lungenparenchym darstellen, und die Dar- 

30 stellung der Lungenparenchympixel in Falschfarben wird die 

Voraussetzung fiir die Befundung von Lungenparenchym geschaf- 
fen. Die Diagnose von PE kann durch die Einbeziehung des Lun- 
genparenchyms erheblich unterstiitzt werden, da Gefafiver- 
schliisse eine Minderdurchblutung bzw. ein Durchblutungsaus- 

3 5 fall im nachgeschalteten Gewebe bewirken und damit eine Ab- 
senkung der Kontrastmittelanreicherung in den entsprechenden 
Bereichen des Lungenparenchyms nach sich Ziehen. Es wird also 
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. deutlich, dass infolge der Falschf arbendarstellung, die im 
Folgenden auch als Farbcodierung bezeichnet wird, die Auswir- 
kungen des thromboembolischen Geschehens direkt sichtbar ge- 
macht werden. 

5 

Im Gegensatz zu der Erf indung sind Veranderungen im Lungenpa- 
renchym in der in der Computertomographie iiblichen Grauwert- 
darstellung aus verschiedenen Grunden nur schwer erkennbar: 

10 Die Verringerung der Kontrastmittelanreicherung wirkt sich 

unter Umstanden in nur geringen Anderungen der HU-Werte aus. 
Dariiber hinaus ist das Lungenparenchym von Gefafien und Bron- 
chien durchzogen. Dabei konnen die GefaSe zu anderen Lungen- 
bereichen gehoren und Kontrastmittel enthalten, wenn ihr Ab- 

15 gang vor dem Gef aSverschluss liegt. Ein Bereich, der durch 

einen Thrombus in dem ihn versorgenden Gefafi homogen minder- 
durchblutet ist, kann deshalb im Bild den Eindruck einer he- 
terogenen Verteilung der HU-Werte vermitteln, mit der Folge, 
dass dieser Bereich durch den Beobachter nicht als mindert 

20 durchblutet wahrgenommen wird. 

Entsprechende Probleme kSnnen im Falle der Erfindung nicht 
auftreten, da diejenigen Bildbereiche^ die kein Lungenparen- 
chym darstellen, also insbesondere die groSen GefaSe in der 
1^5 Lunge, in den gewohnten Grauwerten dargestellt werden. 

Durch die Darstellung in Falschfarben treten auffallig anato- 
mische Bereiche klar hervor. Zum Beispiel sind haufig einzel- 
ne Segmente der Lunge betroffen, die in der Falschf arbendar- 

30 stellung klar sichtbar werden. Da die groSen GefaSe in den 

originalen Grautonen innerhalb des gefarbten Lungenparenchyms 
dargestellt werden, konnen sie vom Radiologen als anatomische 
Landmarken zur Orientierung verwendet werden, was insbesonde- 
re dann von Vorteil ist, wenn die Befundung wie eingangs er- 

35 wahnt anhand von Bilderstapeln erfolgt. 
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. Eine bevorzugte Variante der Erfindung sieht vor, dass als 
erster Verf ahrensschritt durch die Anwendung eines Konturfin- 
dungsalgorithmus in dem CT-Bild die Lungenpixel, d.h. dieje- 
nigen Pixel, die sich innerhalb der zu findenden Kontur be- 
5 finden, von den restlichen Bildinhalten getrennt werden und 
alle folgenden, die Lunge betreffenden Verf ahrensschritte nur 
auf die so identif izierten Lungenpixel angewandt werden. Zur 
Ermittlung der Lungenparenchympixel werden in vorteilhaf ter 
Weise Bronchien und Gefafie arihand ihrer HU-Werte identif i- 

10 ziert und die entsprechenden Pixel aus den Lungenpixeln ent- 
fernt. Dies ist ohne weiteres moglich, da Bronchien HU-Werte 
aufweisen, die unterhalb derer des Lungenparenchyms liegen, 
wahrend die HU-Werte der GefaSe oberhalb der HU-Werte des 
Lungenparenchyms liegen. D.H., wenn die Bronchien und GefaSe 

15 darstellenden Pixel aus den Lungenpixeln entfernt sind, blei- 
ben die Lungenparenchympixel zuriick. Im Folgenden wird die 
Anwendung eines Konturf indungsalgorithmus und Trennung der 
Pixel, die sich innerhalb der zu findenden Kontur befinden, 
von den restlichen Bildinhalten auch als Segmentierung be- 

20 zeichnet- 

Dcunit keine zu grofien L6cher im Bild entstehen, erfolgt die 
Segmentierung gemafi einer Variante der Erfindung so, dass nur 
ein einen vorgebbaren maximalen prozentualen Anteil der Pixel 
,'^^5 nicht uberstei gender Anteil von Pixeln entfernt wird, wobei 
die entfernten Pixel als ungiiltige Pixel eingestuft werden. 

Die Glattungsoperation erfolgt jeweils innerhalb einer wahl- 
baren Filtermaske, d.h. es wird eine Anzahl von Pixeln ge- 

30 wahlt wird, die jeweils in die Glattungsoperation einbezogen 
wird. Dem Pixel im Mittelpunkt der Filtermaske wird der Mit- 
telwert aller in der Filtermaske enthaltenen Pixel zugeord- 
net. Ungiiltige Pixel tragen nicht zum Mittelwert bei , Der er- 
forderliche Mindestanteil von giiltigen Pixeln an den in der 

35 Filtermaske enthaltenen Pixeln kann eingestellt werden. Wird 
der Mindestanteil nicht erreicht, dann wird das Pixel im Mit- 
telpunkt der Filtermaske auf ungultig gesetzt. Damit kann be- 



2001P02362DE 



5 

. stimmt werden, ob die der Filtermaske enthaltenen Locher 
schrumpfen, indem ihre Pixel durch den Mittelwert der sie um- 
gebenden Pixel ersetzt werden, Oder ob sie auf Kosten der 
giiltigen sie umgebenden Pixel wachsen- 

5 

Eine weitere Variante der Erfindung sieht vor, dass der Lun- 
genparenchym darstellende Bildbereich, d.h. die Lungenparen- 
chympixel, einer Glattungsoperation, vorzugsweise einer ins- 
besondere dreidimensionalen gleitenden Mittelung^ unterzogen 
10 wird, was die Befundung erleichtert- 

Eine weitere bevorzugte Ausf uhrungsf orm der Erfindung sieht 
^ vor, dass die in Falschfarben darges tell ten Lungenparenchym- 

pixel dem in Grauwerten dargestellten Originalbild iiberlagert 
15 werden, wobei als ungiiltig eingestufte Pixel durch die ent- 
sprechenden, in Grauwerten dargestellte Pixel des Original- 
bildes ersetzt werden. GroSe GefajSe werden dadurch in ihren 
originalen Grautonen innerhalb der segmentierten und gefarb- 
' ten Lunge dargestellt. 
20 

Urn die Darstellung den jeweiligen Erf ordernissen anpassen zu 
konnen, besteht die Moglichkeit, die in Falschfarben darge- 
stellten Lungenparenchympixel und die in den Grauwerten des 
Originalbildes dargestellten iibrigen Bildbereiche unabhangig 
^5 voneinander Fensterungsoperationen zu unterziehen, Dabei 
sieht eine Variante der Erfindung fiir den Fall, dass eine 
Fensterung der in Falschfarben dargestellten Lungenparenchym- 
pixel erfolgt, vor, dass die Fensterwerte fiir die in Falsch- 
farben dargestellten Lungenparenchympixel aus dem Histogramm 
30 der Lungenpixel ermittelt werden, vorzugsweise indem der 
Schwerpunkt der Histogrammverteilung als Zentral-Wert der 
Fensterung verwendet wird und die Breite des Fensters auf ei- 
nen Wert von ca. 100 HU eingestellt wird. 



35 



Die Anwendung des erf indungsgemaJSen Verfahrens ist nicht auf 
einzelne Originalbilder beschrankt, sondern wird gemaS einer 
bevorzugten Aus fuhrungs form der Erfindung auf die Original- 
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bilder von mehrere Originalbilder umfassende -Volumendaten 
angewendet. wobei auf Basis solcher Volumendaten Bilderstapel 
rekonstruiert werden konnen oder multiplanare Rekonstruktio- 
nen (MPR) vorgenoinmen werden konnen. 

5 

Eine Ausf iihrungsf orm der Erfindung sieht vor, dass zwei ein- 
ander entsprechende Originalbilder angefertigt werden, von 
denen eines ohne und eines mit Kontrastmittel gewonnen wird, 
beispielsweise indem das eine Bild vor und das andere nach 
10 der Gabe von Kontrastmittel aufgenommen wird, und bei dam 

beide Originalbilder nach einem der vorstehend beschriebenen 
Verfahren bearbeitet werden, wobei die beiden erhaltenen be- 
arbeiteten Bilder voneinander subtrahiert werden . 

15 Auch auf diese Weise ist eine gegeniiber den herkommlichen 

Grauwertdarstellungen verbesserte Befundung gewahrleistet . • 

Die Farbdazstellung kann innerhalb des graphischen Userinter- 
faces (Bedienoberf lache) eines CT-Gerats oder einer zur Bear- 
20 iDeitung von CT-Bildern dienenden Workstation als zusatzlicher 
(Soft-) Button (Bedienelement) integriert werden. Bei Betrach- 
ten eines CT-Bildes in beliebiger Orientierung, z.B. nach ei- 
ner Multi Planaren Rekonstruktion (MPR) , kann zwischen der 
Grauwert-Darstellung und der farbigen Lungenparenchym- 
^25 Darstellung. umgeschaltet werden. 

^ Die Erfindung wird nachfolgend an Hand eines in den beige- 
fiigten schematischen Zeichnungen dargestellten Ausftihrungs- 
beispiels erlautert. Es zeigen: 



30 



Fig. 1 ein erf indungsgemaSes CT-Gerat, welches zur Durch- 
fuhrung des erf indungsgemaSen Verfahrens dient. 



35 



Fig. 2 ein die Ermittlung des Startpunktes fur den Kontur- 
f indungsalgorithmus veranschaulichendes Schaubild, 
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Fig. 3 ein den Konturf indungsalgorithinus veranschaulichen- 
des Schaubild, 

Fig. 4 ein die Erosion von auf der gefund'enen Kontur lie- 
5 genden Pixeln veranschaulichendes TSchaubild, wobei 

a) ..den Zustand vor der Erosion und b) den Zustand 
nach der Erosion von 5 lagen von Pixeln darstellt. 



Fig. 5 den Einfluss der der Gef afibeseitigung zugrundelie- 
10 genden Parameter, wobei alle Bilder mit einem Zent- 

ral-Wert von 800 HU und* einer Breite von 300 HU ge- 
fenstert sind und wobei 
a) das Originalbild 

15 
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20 

Fig. 6 den Einfluss des Durchmessers des Filterkems bei 

der Glattiingsoperation wobei der Durchmesser fiir a) 
7 Pixel entsprechend 4.9 mm, fiir b) 9 Pixel entspre- 
chend 6.4 mm und fur c) 11 Pixel entsprechend 7.8 mm 
25 betragt. 

Fig. 7 den Einfluss des Gefafifaktors fur einen Gefafifaktor 
von a) 5 %, b) 20%, c) 35% und d) 50% 

30 Fig. 8 das Userinterf ace des CT-Gerats gemaE Fig. 1, 

Fig. 9 ein axiales Bild der Lunge eines Patienten ohne Be- 
fund in Falschf arbendarstellung, und 

35 Fig. 10 ein axiales Bild der Lunge eines an PE leidenden Pa- 
tienten in Falschf arbendarstellung. 
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• In der Fig. 1 ist ein R6ntgen-CT-Gerat dargestellt, dessen 
insgesamt mit 1 bezeichnete Gantry eine Messoffnung 2 auf- 
weist, welche von einem Drehring 3 iimgeben ist, auf dem ein 
Rontgenstrahler 4 und ein Detektorsystem mit einem in an sich 
5 bekannter Weise bogenformig ausgebildeten, \im eine vorzugs- 
weise parallel zu einer Systemachse Z durch den Fokus des 
Rontgenstrahler s 4 verlaufende Achse gekriimmten, insgesamt 
mit 5 bezeichneten Detektor angebracht sind. Der Detektor 5 
weist mehrere Zeilen 5i bis 5n von Detektorelementen auf, von 
10 denen jede eine Reihe von Detektorelementen auf weist. Von dem 
Rontgenstrahler 4 geht ein strichpunktiert angedeutetes pyra- 
midenf ormiges Rontgenstrahlenbiindel RS aus, dass auf den De- 
tektor 5 auftrifft. 

15 Die Gantry 1 mit dem Rontgenstrahler 4 und dem Detektor 5 ei- 
nerseits und zumindest die zur Aufnahme eines Untersuchungs- 
objektes, z.B. eines Patienten P, vorgesehene Lagerungsplatte 
7 einer insgesamt mit 8 bezeichneten Lagerungsvorrichtung an- 
dererseits sind mittels eines nicht dargestellten motorischen 

2 0 Antriebs in Richtung der parallel zu der Systemachse Z ver- 

laufenden Langsachse der Lagerungsplatte 7 relativ zueinander 
verstellbar, Im Falle des CT-Gerats gemaS Fig. 1 wird dies 
dadurch erreicht, dass die Lagerungsplatte 7 an dem Sockel 9 
der Lagerungsvorrichtung 8 mittels eines Tragteils 10 in 
/^i^5 Richtung der Systemachse Z der Lagerungsplatte 7, also in 

Richtung des mit z bezeichneten Doppelpf eils , verstellbar an- 
gebracht ist. 

Die Lagerungsplatte 7 ist aus einem Rontgenstrahlung wenig 
30 schwachehden Material, z.B. carbonf aserverstarkter Kunststoff 
(CFK) Oder Holz, hergestellt. 

Zur Anfertigung von computer tomographischen Aufnahmen werden 
die Gantry 1 und die Lagerungsplatte 7 relativ zueinander in 
35 eine Position verfahren, in der sich die Lagerungsplatte 7 
durch die Messoffnung 2 der Gantry 1 erstreckt und der auf 
der Lagerungsplatte 7 liegende Patient P eine solche Position 
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. relativ zu der Gantry 1 einninimt, dass ein zu untersuchender 
Bereich des Patienten P von dem Rontgenstrahlenbundel RS er- 
fasst wird, 

5 Zur Anfertigung von Aufnahmen einer oder mehrerer planaren 
Schichten des Patienten P wird der Drehring 3 mit dem Ront- 
genstrahler 4 und dem Detektor 5 zur Aufnahme einer Vielzahl 
von der Rekonstruktion eines oder mehrerer Schnittbilder ei- 
ner Oder mehrerer planare Schichten des Patienten P dienende 

10 Projektionen unterschiedlicher Richtungen urn die Systemachse 
Z gedreht, wahrend die Gantry 1 und die Lagerungsplatte 7 be 
ziiglich der Richtung der Systemachse Z ihre Positionen rela- 
tiv zueinander beibehalten. Da der Detektor 5 mehrere Zeilen 
von Detektorelementen aufweist, konnen gleichzeitig Projek- 

15 tionen beziiglich einer maximal der Anzahl der Zeilen 5i bis 
5n des Detektors 5 entsprechenden Anzahl von Schichten des 
Patienten P aufgenommen warden. Die den Projektionen entspre 
chenden Messwerte, die von einer oder von mehreren Zeilen 5i 
bis 5n des Detektors 5 stammen, sind einem Rechner 11 zuge- 

20 fiihrt, der aus diesen Messwerten in an sich bekannter Weise 
ein Schnittbild, vorzugsweise jedoch einen Stapel Schnittbil 
dem berechnet, die auf einer Anzeigeeinrichtung, z,B. einem 
Monitor 12 dargestellt werden konnen, wobei es fur den fall 
der Erzeugung eines Stapels von Schnittbildern moglich ist, 

4- 

^'1^5 durch den Stapel zu ^blattern*' . 

An den Rechner 11 sind eine Tastatur 13 , eine Mouse 14 
und/oder weitere nicht dargestellte Eingabeeinrichtiangen zur 
Bedienung des CT-Gerats angeschlossen. 

30 

Um die jeweils notwendige Rontgendosis einstellen zu konnen, 
wird die den Rontgenstrahler 4 mit den zu dessen Betrieb er- 
forderlichen Spannungen und Stromen versorgende Generatorein- 
richtung 17 ebenfalls von dem Rechner 11 gesteuert, was iiber 
35 eine entsprechende Leitung angedeutet ist. 
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Zur Durchfuhrung eines sogenannten Spiralscans wird bei akti- 
vem Rontgenstrahler 4 unter kontinuierlicher Rotation des 
Drehrings 3 die Lagerungsplatte 7 geradlinig in Richtung der 
Systemachse Z, also in z -Richtung, verschoben, so dass die 
aufgenornmenen Projektionen sich nicht auf eine oder mehrere 
planare Schichten, sondern auf eine oder mehrere helixfonnige 
Schichten beziehen. Mit Hilfe an sich bekannter Verfahren zur 
Spiralinterpolation ermittelt der Rechner 11 aus den im Zuge 
des Spiralscans gewonnenen gemessenen Projektionen errechnete 
Projektionen, die sich auf eine oder mehrere gewiinschte pla- 
nare Schichten beziehen und die Rekonstruktion entsprechender 
Schnittbilder gestatten. Daruber hinaus ist es moglich, auf 
der Basis von Spiralscans dreidimensionale Bilder zu rekon- 
struieren, da im Zuge eines Spiralscans nicht nur eine pla- 
nare Schicht, sondern ein Volumen abgetastet wird. 

Ftir Untersuchungen mittels eines Kontrastmittels, ist ein 
Kontrastmittelinjektor 15 vorgesehen, mittels dessen liber ei- 
ne Kaniile 16 dem Patienten P ein Kontrastmittel zuftihrbar 
ist. Im Falle des beschriebenen Ausftihrungsbeispiels wird der 
Kontrastmittelinjektor 15, was durch eine entsprechende Lei- 
tung angedeutet ist, von dem Rechner 11 gesteuert, und zwar 
sowohl hinsichtlich der Mehge des Kontrastmittels , die dem 
Patienten P pro Zeiteinheit zugefiihrt wird, als auch hin- 
sichtlich des Beginns und des Endes der Zufuhr von Kontrast- 
mittel . 

Zweckmafiigerweise wird so vorgegangen, dass ziinachst ein Re- 
ferenzscan ohne Aktivierung des Kontrastmittelinjektors 15 
vorgenommen wird, der diejenige Region des Korpers des Pati- 
enten P enthalt, die Gegenstand der Untersuchung ist. 

Die eigentliche Untersuchung wird nach Aktivierung des Kon- 
trastmittelinjektors 15 erst dann gestartet, wenn in dem je- 
weils zu untersuchenden Korperbereich tatsachlich bereits ei- 
ne ausreichende Konzentration des Kontrastmittels vorliegt. 
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im Falle einer Untersuchiing der Lunge im Zusammenhang mit PE 
z.B- nach 30 s (Sekiinden) . 

Das erf indungsgemaSe Verfahren zur Bildbearbeitung der in 
Form von Stapeln von Bildern oder einzelnen Bildem voriie- 
genden CT-Bilder, bei denen es sich Tjim axiale Bilder oder 
durch MPR gewonnene sagittale oder koronale Bilder handeln 
kann, wird von dem Rechner 11 durchgefiihrt und wird im Fol- 
genden naher beschrieben. 

Die Bildbearbeitung erfolgt in 5 Schritten, namlich Segmen- 
tierung, Gef afibeseitigung, Glattung, Farbcodierung und Bild- 
iiber lagerung . 

Zur Segmentierung der Lunge wird ein schwellwertbasierter Al- 
gorithmus verwendet . ^ Der Benutzer kann einen Schwellwert, 
ublicherweise -300 HU vorgeben. Um hohe Funktionssicherheit 
zu erreichen wurde der Algorithmus den speziellen beim Lun- 
genparenchym vorliegenden Bedingungen angepasst. 

Vom Benutzer wird ein primarer Anfangspunkt in beiden Lungen 
gesetzt. Ausgehend von diesem werden mehrere, z.B. sieben, 
sekundare Anf angspunkte ermittelt, wobei drei in y-Richtung 
oberhalb und drei y-Richtxmg unterhalb des primarer Anfangs- 
punkts liegen (Fig. 2) . 

Der Abstand in y-Richtung zwischen den Anf angspunkten betragt 
beispielsweise jeweils 5 Pixel; die x-Koordinaten entsprechen 
denen des primaren Anf ang spunk tes . Aus den Anf angspunkten 
werden in genau horizontaler Richtung, d.h. x-Richtung, sie- 
ben mogliche Startpunkte ermittelt. In jeder Lunge erfolgt 
die Suche nach einem moglichen Startpunkt x-Richtung ausge- 
hend von dem jeweiligen Anfangspunkt in Richtung auf den 
Bildrand. Als Startpunkt wird ein Punkt definiert, dessen CT- 
Wert hochstens gleich dem obengenannten Schwellwert ist und 
der von einer Anzahl von n Pixeln, z.B. n=5 Pixeln, gefolgt 
ist, die oberhalb des Schwellwertes liegen: Der potentielle 
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Startpunkt mit der groSten Entfernung in x-Richtung von dem 
zugehorigen Anfangspunkt wird als der effektiver Startpunkt 
des Konturf indung herangezogen. Dies heiSt, dass jeweils das 
auSerste Pixel als effektiver Startpunkt gewahlt wird, das 
5 noch innerhalb des CT-Wertebereichs fiir Lungenparenchym 
liegt - 

In Schichten mit einer hohen Dichte von mit Kontrastmi tteln 
gefiillten GefaSen wiirde der Standard-Algorithmus mit n=l ge- 

10 legentlich versagen, weil alle Suchpfade an Gefafiwanden und 
nicht am Brustfell (Pleura) enden. Die Annahme eines mittle- 
ren Gef aSdurchmessers von n Pixeln, z.B. n=5, verringert die 
^ Empf indlichkeit des Algorithmus . Der Startpunkt mit der grolS- 
ten Entfernung in x-Richtung von dem zugehorigen Anfangspunkt 

15 wird als der Startpunkt des Konturf indung herangezogen. 

Auch hinsichtlich der eigentlichen Konturf indung der Lunge 
wurde der Standard -Algorithmus angepasst (Fig. 3) . Die Suche 
nach Konturenpunkten geht entgegen dem Uhrzeigersinn vor 
20 sich, ausgehend von dem wie oben beschrieben gefundenen ef- 
fektiven Startpunkt. Der Algorithmus betrachtet immer die 
ersten drei Nachbarpunkte in Suchrichtung und bestimmt zuerst 
das Pixel mit einem Wert unterhalb des Schwellwertes als den 
nachsten Konturenpunkt . Wenn der erste der drei Nachbarpunkte 
^5 als Konturenpunkt detektiert wird, wird die Suchrichtung auf 
minus 90^ gegeniiber der ur sprung lichen Suchrichtung geandert. 
Wenn keines der drei Pixel das Kriterium erfullt, wird die 
Suchrichtung auf plus 90^ gegeniiber der urspriinglichen Such- 
richtung geandert. In alien anderen Fallen wird die ursprling- 
30 liche Suchrichtung unverandert beibehalten. Dem Algorithmus 
ist es gestattet die Suchrichtung umzukehren und somit auf 
seiner eigenen Spur zuriickzugehen . Wenn die Anzahl der Itera- 
tionen einen vorbestimmten Wert iibersteigt, wird die Kontu- 
rensuche abgebrochen . 



35 



Auf Basis der gewonnenen (extrahierten) Kontur wird eine bi- 
nare Segmentierungsmaske erzeugt. Um da^u Pixel zu beseiti- 
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gen, die zum Brustfell gehoren, wird die Segmentierungsmaske 
in fiinf Lagen erodiert, wobei die Erosionen das Prinzip der 
„vier verbundenen Nachbarn'* als Strukturierungselement ver- 
wenden ^ (Fig. 4) . 

5 

Im Falle der Bearbeitung von Stapeln von Bildern, wird die 
Mittelposition jeder Lunge aus der Segmentierungsmaske abge- 
leitet und als primarer Anfangspunkt der Segmentierung des 
nachsten Bildes verwendet . Die Segmentierung schreitet dann 
10 automatisch so lange voran, bis der gesamte Stapel bearbeitet 
ist . 

iU^ Um die Daten fur die folgende Glattungsoperation vorzuberei- 
ten, werden im Zuge der Gef aiSbeseitigung groBere GefaSstruk- 
15 turen und Luftwege durch HU-Wert-Auswahl aus dem Bild ent- 

fernt. Bin unterer Schwellwert HUb und ein oberer Schwellwert 
HUv werden vorgegeben; Pixel unterhalb von HUb werden als 
Luftwege, d.h. Bronchien/ identif iziert , Pixel oberhalb HUv 
als GefaSe. 

20 

Um einen optimalen Bildeindruck zu erhalten, muss eine Balan- 
ce zwischen den beiden Zielen, moglichst alle GefaSe zu ent- 
fernen und moglichst viel Lungenpixel in dem Bild zu behal- 
ten, gefunden werden. Dabei ist der optimale Wert fiir HUv von 
j^25 Patient zu Patient verschieden und kann sich sogar innerhalb 
ein und desselben Patienten andern. Eine Definition von HUv 
als Prozentsatz ist deshalb allgemeingtil tiger , weshalb eine 
Kombination einer schwellwert- und prozentwertbasierten Vor- 
gehensweise Anwendung findet. Untersuchungen haben gezeigt, 
30 dass die maximale Anzahl der entfernten Lungenpixel zweckma- 
Eigerweise auf 28% aller Lungenpixel beschrankt wird, wobei 
HUv so berechnet wird, dass die 28%-Grenze eingehalten wird, 
wahrend fur HUb ist ein fester Wert vori -990 HU zweckmafiig 
ist (Fig. 5) . 



35 



Das segmentierte Bild wird einer Glattungsoperation unterzo- 
gen und dabei im Zuge einer adaptiven Fil.terung durch lineare 
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Interpolation ref ormatiert , um isotrope Pixelabstanden zu er- 
halten. Im Einzelnen kommt eine adaptive gleitende Mittel- 
wertf ilterung unter Verwendung eines isotropen (im 2D-Falle 
kreisf ormigen, im 3D-Falle spharischen) Filterkerns (Kernel) 
zur Anwendung. Wenn wie im Falle des beschriebenen Ausfiih- 
rungsbeispiels keine Lungenparenchymstrukturen unterhalb ei- 
ner GroSe von 5 mm bewertet werden sollen, kommt bevorzugt 
ein Filterkern (Kernel) mit einem Durchmesser von 5 mm ziim 
Einsatz (Fig. 6) . In iiblichen CT-Bildern entspricht dies sie- 
ben Pixeln, 

Bei Anwendung des Filters wird jeweils das mittlere Pixel des 
aktuellen Kernels durch den Mittelwert aller Pixel des jewei- 
ligen Kernels ersetzt. Pixel, die bei vorangegangenen Opera- 
tionen entfernt wurden ( Segment ierung. Erosion, Gef aSbeseiti- 
gung) werden als ungtiltig definiert und tragen zu der Mittel- 
wertbildung nicht bei. 

Durch den Benutzer wird der minimale Anteil der giiltigen Pi- 
xel definiert, die in einem Kernel vorhanden sein miissen, um 
ein giiltiges Ergebnis der Mittelwertbildung produzieren zu 
konnen. Dieser minimale Anteil wird als Prozentwert vorgege- 
ben, der im Folgenden als GefaSfaktor bezeichnet werden wird. 
Wenn der Anteil unterhalb des Grenzwertes liegt, wird das 
mittlere Pixel als ungiiltig gesetzt. Da in dem letzten Verar- 
beitungsschritt alle uncrCiltigen Pixel durch die entsprechen- 
den Pixel des Originalbildes ersetzt werden, legt der GefaS- 
faktor fest, wie viele Gef aSstrukturen und Luftwege in dem 
bearbeiteten Lungenparenchymbereich erscheinen werden, Vor- 
zugsweise ist ein Gefafifaktor von 28% zu verwenden (Fig. 7) . 

Wenn Stapel von Bildern bearbeitet werden, kommt wie bereits 
erwahnt ein spharischer 3D-Filterkern zur Anwendung, der im 
falle eines Durchmessers von sieben Pixeln demnach sieben 
auf einanderf olgende Bilder des Stapels erf asst. 
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Bei der Implementierung des Filteralgorithmus wird vorteil- 
hafterweise von einem schnellen niomerischen Faltungsalgorith- 
mus Gebrauch gemacht . Er wird im Folgenden fiir ein 2D-Filter 
beschrieben und kann leicht auf 3D erweitert werden: 

5 

Alle Pixel des Bildes, die ungultig sind oder auSerhalb der 
gefundenen Kontur liegen werden auf 0 gesetzt. Die Bildmatrix 
und ihre binare Maske werden dann getrennt gefaltet. Die Fal- 
tung der binaren Maske ergibt die Anzahl der gultigen Pixel, 

10 entsprechend der Position des Wertes in der Matrix. Der Ge- 
faSfaktor wird berucksichtigt , indem der entsprechende 
Schwellwert auf die Matrix angewendet wird und alle Werte un- 
f0 terhalb des Schwellwertes auf 0 gesetzt werden. Urn das gefil- 
terte Bild zu erhalten wird das gefaltete Bild durch die ge- 

15 faltete Maske dividiert, wobei die Division Element fiir Ele- 
ment durchgefiihrt wird. Wenn ein Element des gefalteten bina- 
ren Maske 0 ist, wird das Ergebnis auf ungultig gesetzt. 

Urn die Veranschaulichxing der Kbntrastmittelanreicherung im 
20 Lungenparenchym zu erleichtern, wird das resultierende Bild 
in Falschfarben in einer spektralen Farbskala dargestellt, 
wobei eine die Spektralf arben umfassende, also von rot bis 
violett reichende Farbskala zur Farbcodierung zur Anwendung 
kommt. Die zum Zwecke der Falschf arbendarstellung erfolgende 
^25 Farbcodierung wird analog zur Grauwertdarstellung durch eine 
Fensterung gesteuert. Dabei wird zur Ermittlxing der Parameter 
des Fensters heuristisch vorgegangen, d.h. es wird zur auto- 
matischen Bestimmung der Parameter das Histogramm der Lungen- 
parenchympixel analysiert. Dabei wird der Schwerpunkt der 
30 Histogrammverteilung als Zentral-Wert der Fensterxing verwen- 
det und die Breite des Fensters auf einen Wert von ca. 100 HU 
eingestellt. 

Das resultierende Bild des Lungenparenchyms wird defiff Origi- 
35 nalbild liberlagert . Dabei kann eine separate Fensterung der 

in Grauwerten und der in Farbe dargestellten Telle des Bildes 
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erfolgen. Das zugrundeliegende Originalbild ist wesentlich 
fiir die raumliche Orientierung in dem Datensatz- 

Das in Fig. 8 dargestellte Userinterf ace des CT-Gerats gemaS 
5 Fig. 1, das zur Steuerung des erf indungsgemaSen Verfahrens 
dient, umfasst neben einem der Bilddarstellung dienenden Be- 
reich einen Bereich mit Bedienelementen. 

Im einzelnen sind Bedienelemente fiir die Segmentierung vorge- 
10 sehen, wobei mit dem Bedienelement 18 der Schwellwert fiir die 

Segmentierung, mit dem Bedienelement 19 die Anzahl der (pri- 

maren und sekundaren) Anf angspunkte, mit dem Bedienelement 2 0 
^ der Wert fur HUb und mit dem Bedienelement 21 der GefaSfaktor 

eingestellt werden kann. 

15 

Aufierdem ist ein Bereich mit Bedienelementen fiir die Glat- 
tungsoperation vorgesehen, wobei mit dem Bedienelement 22 der 
Filterkern samt Durchmesser des Filterkerns und mit dem Be- 
dienelement 23 im falle der Bearbeitung eines Stapels von 
2 0 Bildern die Anzahl der von dem Filterkern erfassten Bilder 
gewahlt wird. 

Weiter sind zwei Bereich 24 und 25 vorgesehen, die Bedienele- 
mente zur Einstellung der Parameter der Fensterung der 
5 schwarzweiSen und farbigen Bildanteile enthalten, 

Mit Bedienelementen 27 und 2 8 besteht die Moglichkeit zwi- 
schen SchwarzweiS- und Falschf arbendarstellung \imzuschalten. 

30 SchlieSlich ist ein Bereich 29 vorhanden, der Bedienelemente 
enthalt, die dazu dienen zwei nach dem der vorstehend be- 
schriebenen Verfahren bearbeitete Bilder voneinander zu sub- 
trahieren, von denen eines ohne und eines mit Kontrastmittel 
gewonnen wurde. 

35 

Der erf indungsgemaSen Bildbearbeitung wurden klinische Daten- 
satze unterzogen, die bei Brustuntersuchungen aufgenoramen 
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wurden, die wegen klinischen Verdachts auf PE durchgefiihrt 
warden. Dabei wurde nichtionisches Kontrastmittel (Ultravist 
370, Schering, Berlin, Germany) intravenos mit einem Kon- 
trastmittelinjektor (CT 9000 Digital Injection System, Lie- 
5 bel-Flarsheim, Cincinnatti/ OH) verabreicht . Dabei wurden ins- 
gesamt 120 cm^ Kontrastmittel bei einer Flussrate von 3 cm"^/s 
verabreicht, gefolgt von einer Kochsalzlosung (insgesamt 
30 cm*^ bei 3 cm^/s) . Die Startverzogerung zwischen dem Beginn 
der Kontrastmittelgabe und dem Beginn der Aufnahme betrug bei 

10 alien Patienten 30 Sekunden. Die Aufnahme der Bilder erfolgte 
mit einem handelsiiblichen Mehrzeilen-CT-Gerat (Somatom Volume 
Zoom; SIEMENS AG, Forchheim, Germany) . Die Auf nahmeparameter 
>^ waren 140 kV und 100 mAs unter Verwendung einer engen Kolli- 
mierung von 4x1 mm bei einem Pitch von 1.75. Daher konnte 

15 die gesamte Brust wahrend einer einzigen Atempause von 21 s 
unter sucht werden. Es waren keine weiteren Unter suchungen in 
spateren Atempausen erf order lich. Die so gewonnenen Ausgangs- 
bilder wurden retrospektiv fiir eine Schichtdicke von 1,25 mm 
und ein Rekonstruktionsinkrement von 0,8 mm rekonstruiert und 

20 zur Bildbearbeitung auf einen PC (Pentixim III, 600MHz) expor- 
tiert. 

Um die Wirkung des erf indungsgemaSen Verfahrens zu demonst- 
rieren, sind Bilder eines Patienten ohne pathologische Veran- 
^^.5 derung (Fig. 9) und eines Patienten mit erwiesener PE (Fig. 

10) jeweils in axialer Orientierung dargestellt, wobei in der 
Praxis die Bereiche hoher Dichte, d.h. hoher Kontrastmittel- 
anreicherung, rot und die Bereiche geringer Dichte, d.h. im 
wesentlichen ohne Kontrastmittelanreicherung, violett darge- 
3 0 stellt sind. Wegen der Schwarzweifi-Darstellung der Figuren 9 
und 10 sind in diesen die Bereiche hoher Dichte durch mit H 
und die Bereiche geringer Dichte durch mit L bezeichneten 
Pfeilen gekennzeichnet . 

35 Fig. 9: Axiales Bild auf der Hohe des linken Atriums nach der 
intravenosen Gabe von Kontrastmittel. Es liegt eine vollstan- 
dige Kontrastmittelanreicherung in alien Lungengef aSen vor. 
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Es sind keine Perfusions-Def izite in der Darstellung des Lun- 
genparenchyms sichtbar. Es sind weder Bereiche besonders ho- 
her noch besonders geringer Dichte erkennbar, weshalb ent- 
sprechende Pfeile fehlen. 

Fig. 10: Mehrere Arterien weisen Verschliisse durch Thromben 
auf . Entsprechend treten Bereiche erhohter Dichte durch ge- 
stautes Kontrastmittel und Bereiche verringerter Dichte in- 
folge Minderdurchblutung auf, die durch entsprechende Pfeile 
gekennzeichnet sind . 

Es wird also deutlich, dass mittels des erf indungsgemaSen 
Verfahrens eine verbesserte Befundung moglich ist. Beid Enden 
der Farbskala, d.h. rot und violett, konnen auf eine patholo- 
gische Situation im Lungenparenchym hindeuten. Rote Bereiche 
sind Bereiche erhohter Dichte und konnen beispielsweise durch 
Atelektasen verursacht werden. Violet te Bereiche sind Berei- 
che verringerter Dichte, die moglicherweise durch verringerte 
Durchblutung und damit verringerte Kontrastmittelanreicherung 
verursacht sein konnen. 
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Patentanspriiche 

1. Verfahren zur Bildbearbeitung ausgehend von einem unter 
Verwendung eines Kontrastmittels gewonnenen Computertomogra- 
phie (CT) -Bildes einer Lunge als Originalbild, aufweisend die 
Verf ahrensschritte : 

Ermittlung der Lungenparenchympixel , 

Darstellung der Lungenparenchympixel in Falschf arben, und 

Darstellung der tibrigen Bildbereiche in den Grauwerten des 
Originalbildes . 

2. Verfahren nach Anspruch 1, bei dem als erster Verfahrens- 
schritt durch Anwendung eines Konturf indungsalgorithmus in 
dem Originalbild die Lungenpixel von den rest lichen Bildin- 
halten getrennt werden und alle folgenden die Lunge betref- 
fenden Verf ahrensschritte nur auf die so identif izierten Lun- 
genpixel angewandt werden. 

3- Verfahren nach Anspruch 1 oder 2, bei dem zur Ermittlung 
der Lungenparenchympixel Bronchi en und GefaJSe an Hand ihrer 
HU Wert identif iziert und die entsprechenden Pixel aus den 
Lungenpixeln entfernt werden. 

4 . Verfahren nach Anspruch 3 , bei dem ein einen vorgebbaren 
maximalen prozentualen Anteil der Pixel nicht libersteigender 
Anteil von Pixeln entfernt wird. 

5. Verfahren nach Anspruch 4, bei dem die entfernten Pixel 
werden als ungiiltige Pixel eingestuft werden. 

6. Verfahren nach einem der Anspruche 1 bis 5, aufweisend den 
Verf ahrensschritt , dass der Lungenparenchym darstellende 
Bildbereich einer Glattungsoperation unter zogen wird, wobei 
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ungtiltige Pixel nicht in die Glattungsoperation einbezogen 
werden . 

7. Verfahren nach Anspruch 6, bei dem als Glattungsoperation 
eine gleitende Mittelung vorgenommen wird. 

8 . Verfahren nach Anspruch 6 oder 7 , bei dem eine Anzahl von 
Pixeln gewahlt wird, die jeweils in die Glattungsoperation 
einbezogen wird. 

9. Verfahren nach Anspruch 8, bei dem die Anzahl von jeweils 
in die Glattungsoperation einbezogenen Pixeln ein mittleres 
Pixel umfasst, dem Pixel der Mittelwert aller in die Glat- 
tungsoperation einbezogener Pixel zugeordnet wird. 

10. Verfahren nach Anspruch 9, bei dem ein Mindestanteil von 
gultigen Pixeln an der Anzahl von jeweils in die Glattungs- 
operation einbezogenen Pixeln gewahlt wird, bei dessen Unter- 
schreitung das mittlere Pixel als ungiiltig gesetzt wird. 

11. Verfahren nach einem der Anspruche 1 bis 10, bei dem die 
in Falschfarben darges tell ten Lungenparenchympixel dem in 
Grauwerten dargestellten Originalbild uberlagert werden, wo- 
bei als ungrultig eingestufte Pixel durch die entsprechenden 
in Grauwerten dargestellten Pixel des Originalbilds ersetzt 
werden . 

12. Verfahren 1 bis 11, bei dem die in Falschfarben darge- 
stellten Lungenparenchympixel und die in den Grauwerten des 
Originalbildes dargestellten iibrigen Bildbereiche unabhangig 
voneinander Fensterungsoperationen unterzogen werden. 

13. Verfahren nach einem der Anspruche 1 bis 12, bei dem eine 
Fensterung der in Falschfarben dargestellten Lungenparenchym- 
pixel erfolgt und die Fensterwerte ftir die in Falschfarben 
dargestellten Lungenparenchympixel aus dem Histogramm der 
Lungenpixel ermittelt werden. 
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14. Verfahren nach Anspruch 13, bei dem der Schwerpunkt der 
Histograiran-Verteilung wird als Zentral-Wert der Fensterung 
verwendet wird und die Breite des Fensters auf einen festen 
Wert von ca, 100 HU eingestellt wird. 

15. Verfahren nach einem der Anspriiche 1 bis 14, welches auf 
die Originalbilder von mehrere Originalbilder umf assenden Vo- 
lumendaten angewendet wird. 

16. Verfahren nach Anspruch 15, bei dem auf Basis von mehrere 
Originalbilder umfassenden Volumendaten eine multi-planare 
Rekonstruktion (MPR) vorgenommen wird. 

17. Verfahren, bei dem zwei einander entsprechende Original- 
bilder angefertigt werden, von denen eines ohne und eines mit 
Kontrastmittel gewonnen wird, und bei dem beide Originalbil- 
der nach einem der Verfahren nach einem der Anspriiche 1 bis 
16 bearbeitet werden, aufweisend den Verf ahrensschritt , dass 
die beiden erhaltenen bearbeiteten Bilder voneinander subtra- 
hiert werden. 

18. CT-Gerat, welches wenigstens ein Bedienelement zur Durch- 
fxihrung eines Verfahrens nach einem der Anspruche 1 bis 17 
aufweist • 

19. CT-Gerat nach Anspruch 18, bei dem das Bedienelement Be- 
standteil eines graphischen Userinterf aces ist. 

20. Workstation zur Bearbeitung von CT-Bilder, welche wenigs- 
tens ein Bedienelement zur Durchfiihrung eines Verfahrens nach 
einem der Anspruche 1 bis 17 aufweist. 

21. Workstation nach Anspruch 20, bei dem das Bedienelement 
Bestandteil eines graphischen Userinterf aces ist. 
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Zusairanen fas sung 

Verfahren zur Bildbearbeitung ausgehend von unter Verwendung 
eines Kontrastmittels auf genoiranenen Computertomographie (CT) - 
Bildes einer Lunge und CT-Gerat zur Durchfuhrung eines sol- 
Chen Verfahrens 

Die Erfindung betrifft ein Verfahren zur Bildbearbeitung aus 
gehend von einem vorzugsweise unter Verwendung eines Kon- 
trastmittels auf genoiranenen Computertomographie (CT) -Bildes ei 
ner Lunge als Originalbild. Dabei sind als die Verfahrens- 
schritte die Ermittlung der Lungenparenchympixel , die Dar- 
stellung der Lungenparenchympixel in Falschf arben, und die 
Darstellung der ubrigen Bildbereiche in den Grauwerten des 
Originalbildes . 
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